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ABSTRACT 



We reevaluate the hadronic part of the electromagnetic vacuum expectation 



■^ value using the standard dispersion integral approach that utilizes the hadronic 

■^ cross section measured in e~^e experiments as input. Previous analyses are based 

^ upon point-by-point trapezoidal integration which has the effect of weighting all 

o. 



(]j inputs equally. We use a technique that weights the experimental inputs by their 



stated uncertainties, includes correlations, and incorporates some refinements. We 
find the hadronic contribution to the fractional change in the electromagnetic cou- 
pling constant at q^ = M|, Aq!(M|), to be 0.02666 ± 0.00075, which leads to a 
value of the electromagnetic coupling constant, a~^{M'^) = 129.08 ± 0.10. This 
value significantly shifts the Standard Model predictions for the effective weak 
mixing angle measured at the Z pole and moderately shifts the predicted Z width. 
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1. Introduction 

At the current time, a large program of precise electroweak measurements is 
being conducted throughout the world. The object of this program is to test 
the electroweak Standard Model by comparing the measured values of a large set 
of electroweak observables with the predictions of the Minimal Standard Model 
(MSM). The Standard Model calculations have been performed to full one-loop 
accuracy and partial two-loop precision by a large community of researchers. In 
all of these calculations, it is necessary to evaluate the one-particle-irreducible 
contributions to the photon self-energy Il^^{q'^) or the related quantity H' (q^) = 
(H-y-ylq ) — ll^-y{0))/q 0.1 the Z mass scale q = MJ. These quantities are usually 
absorbed into the definition of the running electromagnetic coupling a{q'^), 



a(q ) =- 



"i-n;,((?2)' 

where aQ = 1/137.0359895(61) is the electromagnetic fine structure constant. This 
quantity is also represented as the fractional change in the electromagnetic coupling 
constant Act, 

Using analytic techniques and the optical theorem applied to the amplitude for 
s-channel Bhabha scattering, the quantity Aa has been related to the cross section 
for the process e"'"e~ — >■ 7* — )■ all {(Ttot) as follows, 

00 

where Rtot{s) = <ytot{s) I <y ^ij.{s) is the ratio of the total cross section to the (mass- 
less) muon pair cross section at the center-of-mass energy ^/s. It should be noted 
in passing that equation (3) is correct to all orders in tto and relies only upon 



the assumption that the real part of H^^ is much larger than its imaginary part 
(the next-order correction is proportional to Im n^^/|n^^|^ which is approximately 
3xl0~ at (f = M^)- It is straightforward to evaluate equation (3) for the con- 
tinuum leptonic cross sections. In the limit that the scale q^ is much larger than 
the square of the lepton mass m^, the contribution of the continuum leptonic cross 
sections is given by the following expression, 



A«,(./) = |^E 



-- + ln— 
3 rnj 



(4) 



The remaining contributions to Rtot consist of the continuum hadronic cross 
section and the J = l~ resonances and are labelled Rhad- Since the cross sections 
for these final states are not accurately calculable from first principles, experimental 
inputs are used to evaluate the remaining portion of equation (3), 



/2 
'^^~r^ \^had{s). (5) 

s{q^ - s) 
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Equation (5) has been evaluated at the Z boson mass scale several times. 
A complete description was given by Burkhardt, Jegerlehner, Penso, and Verzeg- 
nassi in 1989. The result was updated by Jegerlehner in 1991 (to include 
measurements from the Crystal Ball Collaboration), 

r 0.0288 ± 0.0009, Reference 4 

AahadiMz) = \ (6) 

I 0.0282 ±0.0009, References. 

The authors of Reference 4 perform the integration in three parts: the hadronic 
continuum above ^/s = 1 GeV, the 7t'^7i:~ final states above threshold; and the 
u, (j), ij), and T resonances. The continuum integration is performed by linearly 
interpolating between the data points. The resonance contributions were calculated 
from an analytic expression which results from integrating a Breit-Wigner lineshape 



and depends upon the masses, widths, and leptonic widths of each resonance. The 
uncertainties on each contribution were estimated by techniques which appear to 
be very conservative. The work reported in this paper was undertal^en initially 
to estimate a more accurate error. We have indeed performed what we believe 
to be a more accurate analysis, however, it is our conclusion that the uncertainty 
estimated by the authors of Reference 4 is not overestimated. 

2. The Analysis 

Any attempt to combine the results of many experiments is a perilous un- 
dertaking. Many different techniques and approaches have been used. Not all 
researchers have addressed all possible problems nor are systematic error estimates 
performed in uniform ways or to uniform standards. We therefore adopt some 
the techniques of the Particle Data Group. Older measurements which are con- 
tradicted by newer, more precise work are excluded from the analysis. Parameter 
uncertainties that are extracted from fits with x^ per degree of freedom (dof) larger 
than one are rescaled by the factor \/x^ jdof . 

2.1 The Data 

The approach to the evaluation of equation (5) is driven by the form of the 
data themselves. The total hadronic cross section can be decomposed into four 
pieces: the hadronic continuum above W = ^/s = 1 GeV, the charged two-body 
final states tt+tt" and K'^K~ from their respective thresholds to 2.6 GeV, and 
hadronic resonances (excluding charged two-body final states). Since equation (5) 
is linear in the hadronic cross section, we decompose Aa^ad as follows, 

Aaw(?') = A^ril^^) + AaCf{q') + Aafj-{q^) + AaTM), (7) 

where the four terms on the right-hand side correspond to the four pieces of the 
hadronic cross section. 



The rationale for this decomposition is as follows. The region below W = 
1 GeV is dominated by the p, co, and (p resonances. The electromagnetic form 
factors for the processes e"'"e~ — )■ 7t'^7t~ and e"'"e~ — )■ K'^ K~ are mea- 

sured well from threshold to W^ ~ 2 GeV. Resonances do not account for all of 
the TT+vr" and K^K~ cross section in this region. On the other hand, essentially 

all other two-body and three-body final states are associated with the resonances. 

ri9i 
Measurements of three-pion final states near W = 1 GeV show the non-resonant 

portion to be consistent with zero. Similarly, measurements of various two-body 
final states such as KiKg show small non-resonant cross sections. The cross sec- 
tions for four-pion final states become significant above 1 GeV but are small below 
that energy. The 772 experiment at the ADONE storage ring at Frascati has 
measured the hadronic cross section ratio for three or more hadron final states, Rr , 

' had 

from W = 1.42 GeV to M^ = 3.09 GeV. They have also presented several points 
from 1 GeV to 1.4 GeV that are composed of various multipion cross sections from 
Novosibirsk and Orsay ' ' and are claimed to approximate i?^^^. Measurements 
beginning at M^ = 2.6 GeV by the MARK I,"'^ DASP,"'^ and PLUTO "'^ Collabo- 
rations claim to measure the entire cross section. We therefore conclude that Rhad 
is well approximated below Wi = 2.6 GeV by a sum of the 7t'^7t~ and K'^ K~ 
contributions from threshold to W\ (where they are much smaller than Rj^^^ ; the 
^had measurements from 1 GeV to W\] and the p, uj and (/) resonances where the 
hadronic widths are adjusted to remove the 7t'^7t~ and K'^K~ final states that are 
already included explicitly. Note that the several broad e"'"e~ resonances between 
the (/)(1020) and W = 2 GeV are implicitly contained in the two-body or Rj^^^ 
categories. Since the tt+tt" and K^K~ cross sections are very small at W\, the 
^had ^^^ total continuum Rhad measurements should be continuous at this point. 

At center-of-mass energies larger than W\, many measurements of the hadronic 
continuum and resonances exist. The region of the charm threshold from W = 
3.6 GeV to M^ = 5.0 GeV is complicated and not well measured. The MARK I, 
DASP, and PLUTO Collaborations all observe an enhancement beyond the ex- 
pected threshold shape. The DASP data show three resolved resonances. The 



MARK I and PLUTO data are consistent with the DASP data but do not cleanly 
resolve the resonances. We choose to follow the Particle Data Group and recognize 
the DASP resonances: ^(4040), -0(4160), and ^'(4415). The ip family therefore 
consists of six states. 

Between 5 GeV and 10.4 GeV, the MARK I, DASP, PLUTO, Crystal Ball,"'^ 
LENA,'"^ CLEO,"'^ CUSB,"'^ and DESY-Heidelberg"°^ Collaborations have pub- 
lished Rhad measurements which are are plotted in Figure 1. The error bars include 
only point-to-point uncertainties. The recently published Crystal Ball measure- 
ments have a systematic normalization uncertainty of 5.2%. The other measure- 
ments have normalization uncertainties in the range 6.8-10%. The data are also 

[311 

compared with the recent QCD prediction of Chetyrkin and Kuhn which in- 
cludes quark mass effects. At VU = 5 GeV, the MARK 1 data are consistent with 
other measurements. As W increases, they show a systematic increase in Rhad 
and suggest the presence of a structure near 6.6 GeV. Including the quoted 10% 
normalization uncertainty, the MARK 1 data are larger than the more precise mea- 
surements by approximately two standard deviations. The reader is reminded that 
first generation detectors like MARK 1, DASP, and PLUTO were small acceptance 
devices that necessarily involved large acceptance corrections without the benefit 
of good event structure modelling. After acceptance corrections and a r-lepton 
subtraction, the MARK 1 group observed that two-charged-prong events consti- 
tuted nearly 20% of the hadronic cross section of i? at VU = 7 GeV. This is about 
1.5 times the two-prong rate due to r'^T~ final states and three times the rate that 

[32l r-n m [33l 

is predicted by the JETSET 7.3 Monte Carlo program. While this may not 
be wrong, we choose to exclude data from the first generation experiments when 
more modern results are available. Such data are available above charm threshold. 
Unfortunately, we are constrained to use very old continuum measurements below 
charm threshold. 

The Particle Data Group lists six T family resonances between 9.4 GeV and 
11 GeV. All are included in the resonance contribution. 

6 
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Figure 1. The Rhad measurements of the MARK I (circles), PLUTO 
(horizontal marks), Crystal Ball (diamonds), LENA (squares), CLEO 
(stars), CUSB (X's), and DESY-Heidelberg (vertical crosses) Collaborations 
in the region between W = 5 GeV and W = 9.4 GeV. The error bars include point- 

[311 

to-point uncertainties only. A recent QCD calculation which includes quark mass 
effects is shown as a solid line for as{M^) = 0.125. 

Above b-quark threshold, a number of Rhad measurements have been car- 
ried out by the PEP and PETRA experiments. However at energies above 
W = 34: GeV, Z-7 interference becomes significant. We therefore use only those 
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measurements in the region W < 34 GeV where the correction for electroweak 
interference is less than 1%. 

We expect that Rhad is well described by perturbative QCD in the region 
above b-quark threshold. This implies that the measurement of the strong coupling 
constant, as{M'^), extracted from a fit to the lineshape parameters of the Z boson 
by the LEP Collaborations, 

a^(M|) =0.125 ±0.005, (8) 



provides a precise measurement of Rhad at M^ = Mz- We have verified that the 
value of as{M'^) given in equation (8) is insensitive to the value of the electro- 
magnetic coupling constant a{M'^) used in the fitting procedure (tt(M|) can be 
left as a free parameter with essentially no effect on the extracted value of as). 
To convert as{M'^) into a determination of Rhadi^z)^ we use the third-order MS 



QCD expression. 



[41] 



^ o (3-/9?) 

RQCD{s)=3Y,Q)Pf^Y^ 
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(9) 



where: Qf is the final state fermion charge, (]f = a 1 — Arn^J s is the fermion veloc- 
ity in the e"'"e~ center-of-mass frame {nif is the fermion mass), and the coefficients 



are functions of the number of active flavors N 



/' 



r2 



1.9857 -0.1153iV 



/ 



6.6368 - 1.2002iVf - 0.0052iV? - 1.2395 



'f 



f 






The resulting value of Rhadi^z) is, 



(10) 



i?w(Mz) = 3.818 ±0.006. 



(11) 



The following three sections of this chapter describe the evaluation of: the 



continuum contribution Aa'j^^ , the contributions of the charged two-body final 
states Aal J and Aa^ /^ , and the resonance contribution Aa]^^,. 

had had ' had 

2.2 The Hadronic Continuum 

The authors of Ref. 4 evaluated equation (5) for the continuum contribution 
by performing a trapezoidal integration with measured values of Rhad- Their ap- 
proach has two advantages: it is unbiased by human prejudice about the functional 
form of Rhad(s), and it would automatically account for undiscovered resonances 
which are broad as compared with the spacing of measurements. Unfortunately, 
this technique also has a serious shortcoming: it does not take experimental errors 
into account properly. All data points receive equal weight irrespective of their 
experimental precision. An experiment which publishes a large number of impre- 
cise data points receives more weight than an experiment which publishes fewer 
precisely measured ones. 

We avoid this problem by fitting the data to an appropriate functional form 
Rfit{s; a^,) where a^, are the parameters of the function. In the absence of undiscov- 
ered resonances, Rhad can be described by a continuous function. A x^ fif has the 
virtue that measurements are weighted by their experimental errors. To do this, 
we make the conservative assumption that all normalization uncertainties within 
an appropriate grouping of measurements are 100% correlated. The x function 
therefore has the form. 



x' 



= Yl [^had - Rf^t{s^] A/o) W^^' Rhad " Rfrt{Sj]aT, 



fl2) 



h3 



where i?|j^j is the value of Rhad measured at energy Sj and the inverse elements of 
the weight matrix Wjj are given by the following expression, 

{af (ptp) + af (norm) , i = j 
(Ti (norm) (Tj (norm), ^ 7^ J, same grouping (13) 

0, ^ 7^ i, different grouping 

where cTj (ptp) and (Tj(norm) are the point-to-point (statistical and systematic) and 



normalization uncertainties associated with the i measurement. 

Equation (5) is evaluated by performing a Simpson's Rule integration using 
the function Rfn and the best estimate of the parameters. The parameter un- 
certainties 5ak reflect the point-to-point and normalization uncertainties to some 
extent. Unfortunately, the process of fitting a large number of measurements with 
a function of a smaller number of parameters necessarily involves some loss of in- 
formation. If we add a priori information to the problem by choosing a physically 
motivated fitting function, the information contained in the parameter error matrix 
may be entirely appropriate. To understand this problem better, we evaluate the 
uncertainty on Aa^adiM^) by two techniques. In the first, the parameter uncer- 
tainties are propagated to the calculated value of Aa^adiM'^) using the following 
expression which is valid for any function of the parameters, 



^^i^a,.^,U, = j:'-^^^E„'-^^^. (14, 



where the derivatives are calculated numerically and Ef^i = (Saj^Sai) is the param- 
eter error matrix that is extracted from the fitting procedure. The second error 
estimate is performed by constructing a large ensemble of data sets by shifting the 
measured data points i?|j^j(meas) as follows, 

i?Ld(set j) = i?Ld(meas) + /,f a,(ptp) + /,7-a, (norm) , (15) 

where the factors f^j are Gaussian-distributed random numbers of unit variance. 
The entire fitting and integration procedure is then applied to each member of the 
ensemble. The uncertainty on Aa^adiM^) is determined from the central 68.3% 
of the ensemble distribution. 

The use of a fitting function has the problem that one may introduce bias 
through the choice of parameterization. We attempt to evaluate this effect by 
varying the parameterizations as much as ingenuity and computer time allow. The 
quoted contributions to Aa^adiMz) ^^^ those corresponding to the best fits. Each 
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contribution is assigned a parameterization uncertainty S{Aahad)param based upon 
the spread of results corresponding to reasonable fits. 

The first step in the evaluation of equation (5) for the hadronic continuum 
is to formulate a suitable (piecewise-continuous) parameterization Rf^t{s]a}.). We 
choose to use the perturbative QCD expression given in equation (9) in the region 
M^ > 15 GeV and an empirical parameterization in the region 1 GeV< W < 
15 GeV. In the high energy region, the only free parameter is as{M'^) which is 
evolved to other scales using the prescription given by Marciano. 

In the portions of the low energy region that are measured well, polynomials are 
used to parameterize R^adi^)- To ensure that the function is continuous across 
several points Wa, the polynomials are constructed in Xa = W — Wa and the zeroth 
order terms are excluded, 

n 

1=1 

where a is a label to distinguish different regions. Separate polynomials are used 
to describe the following regions: 1 GeV< ly < 1.9 GeV (labelled region s), 
1.9 GeV< W < 3.6 GeV (labelled region c), and 5.0 GeV< W < 10.4 GeV (labelled 
region b). Although a single, large-order polynomial is adequate to describe the 
data between W = 1 GeV and charm threshold at 3.6 GeV, the data show a 
distinct shape change near W = 1.9 GeV (where the four-pion cross section is 
becoming small). It was possible to obtain better fits by introducing an additional 
polynomial to describe the region from 1 GeV to 1.9 GeV. A comparison of the two 
possible forms is used to assess the parameterization sensitivity of the final result. 

Since there are no measurements of the continuum Rhad in fh^ b-quark and 
c-quark threshold regions, it is necessary to extrapolate the form of Rhad from 
3.6 GeV to 5.0 GeV and from 10.4 GeV to 15 GeV with functions that are physi- 
cally motivated. In the case of the charm threshold region, the DASP Collaboration 
has published (in graphical form) the shape of the continuum that was preferred 
by their resonance fits. The function which characterizes the shape of the thresh- 
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old, fDASPi^)? does not increase as sharply as the free-quark threshold factor 
(3(3 — /3^)/2 but increases more rapidly than the (3^ threshold factor for pointlilce 
scalar particles. To construct the function Rfn, all three possibilities are used for 
the c-quark threshold and the two extreme possibilities are used for the b-quark 
threshold, 

fc{W) = { fDASp(W) hiW) = \'\ ^ ^' (17) 

where the c- and b-quark masses are taken to be the D and B meson masses, 
respectively. The actual size of the charm-associated step in Rhad, ^Rc is left as 
a free parameter. The size of the bottom-associated step in Rhad is constrained to 
be the difference between the value of the fit function at W = 10.4 GeV and the 
value of the QCD portion at 1^ = 15 GeV, AR^ = RqcdC^^) - i?/,i(10.4). 

The actual form of the fitting function is given by the following expression, 

( Ro + P^^{W -1.0), l<iy<1.9 

i?/,f(1.9) + P^^(M/- 1.9), 1.9<W<3.6 

Rfu{3.6) + ARJ,{W), 3.6<M/<5.0 

Rfu{5.0) + P^^{W -5.0), 5.0 < ly < 10.4 

RfuilOA) + ARhhiW), 10.4 <W< 15.0 

{Rqcd(W), 15 <W 

where Rq, the value of Rhad at W^ = 1 GeV, was a free parameter and the order of 
the polynomials varied from 1 to 7. The number of free parameters varied from 7 to 
30. The fit quality did not improve substantially when the number of parameters 
exceeded 10. The weight matrix was constructed from equation (13) assuming that 
normalization uncertainties were completely correlated in four groups: the 20% 
uncertainties of the lowest energy measurements (1.0 GeV< W < 3.09 GeV), 

the 15-20% uncertainties of the MARK 1, DASP, and PLUTO measurements 
(2.6 GeV< W < 4.9 GeV), the 5-10% uncertainties of the measurements 
between charm and bottom thresholds, and the 1.7-7.0% uncertainties of the PEP 
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RfuiW) 



and PETRA experiments above bottom threshold. The data are corrected 



for electroweak interference before the fitting procedure is applied. The data and 
the result of a typical fit are shown in Fig. 2. The error bars include the point-to- 
point and the normalization uncertainties. The fit quality was excellent (x^/dof = 
104.7/105). 
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Figure 2. The continuum Rhad measurements including normalization uncertainties. A typ- 
ical fit to equation (18) is shown as the solid curve. The dashed curve corresponds to the 
equal-weighting test described in the text. 

The various hypotheses for Rfn are used to evaluate the integral in equation (5) 
from sq = 1 GeV^ to oo = 10^ GeV^. Although the singularity in the integrand is 
formally well controlled, digital computers are famous for their inability to under- 
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stand formalities. We have therefore recast equation (5) into a form which is more 
suitable for electronic evaluation, 



g -A 

2 aoq'^ \ RMq'^) . \q^-so] f Rf^t{s) - RfrM^) 

^ahad[q ) = ^^{ -^-^ In - ds- ' 



StT I fl2 g^ j g(g — q2 



dR 



fit 



ds 
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-So' 
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50 
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+ A" 
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-A. 



oo 



(19) 



In —z — ds 



Rf^t{s) - Rftt{q 



'2\ 



-,2 rf^ — 1\ I at c — n'2\ 



S[S 



where we have assumed that Rf^t is well approximated by a linear expansion over 
the interval q^ — A < s < cp' + A. (in practice, we use A = 0.5 GeV^). The evaluation 
of equation (19) requires that Rqcd(s) be extrapolated through t-quark threshold. 
For this purpose, the top quark mass is assumed to be 175 GeV. 

The contribution of the hadronic continuum to Att/jad(Mj) is found to be fairly 
insensitive to the form of Rj^t- The central value of Aaf^^^lM'^) corresponds to the 
best estimate of the parameters of the function which uses: the DASP shape for the 
c-quark-threshold, the free-quark shape for the b-quark-threshold, and the values 
(2,3,3) for {Ni,,Nc,Ns). The maximum deviation from this value occurs when the (3^ 
function is used for the c-quark-threshold and a sixth-order polynomial is used to 
parameterize the entire region W < 3.6 GeV. The size of the maximum deviation 
is taken as an estimate of the parameterization uncertainty. The experimental 
uncertainty given by equation (14) is found to be a smaller by a factor of three 
than the estimate derived from the ensemble of fluctuated data sets. Since our 
fitting function has no physical motivation whatsoever in the low W region, we 
choose the larger estimate as the more accurate. The resulting contribution to 

AahadiMz) is 

^«w (^1) = 0.021428 ± 0.000724(exp) ± 0.000150(param). (20) 

The experimental uncertainty given in equation (20) dominates the uncertainty 
on Aa)iad{M'^). We note that this uncertainty is itself dominated by the 0.000655 
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contribution of the normalization uncertainties of the Rhad measurements below 
charm threshold (particularly in the region 2.6-3.6 GeV). Any further improve- 
ment in the uncertainty on AahadiMz) requires that improved measurements be 
performed in this region. 

The central value of this result is somewhat smaller than the one given in 
Ref. 4 for four reasons. The first is that we've defined the continuum contribution 
to exclude the charm-related enhancement near 4 GeV. The charm-threshold re- 
lated enhancement is incorporated by the inclusion of the '0(4040), '0(4160), and 
'0(4415) resonances in the resonance contribution. The second difference is that 
our technique weights input data by their uncertainties and accounts for the large 
correlated uncertainties between the measurements within a measurement group. 
The third difference is that we have replaced the 44 MARK I measurements of the 
continuum between W = 4.9 GeV and W = 7.6 GeV with the more recent Crystal 
Ball data. The fourth difference is that we use the LEP measurement of as{M'^) 
to constrain Rhad in f h^ high energy region to a somewhat smaller value than the 
one preferred by the PEP/PETRA experiments alone. To verify that these differ- 
ences are indeed source of the discrepancy, we have repeated the analysis with all 
points weighted equally (all points have constant fractional uncertainties and are 
assumed to be uncorrelated). To ensure that the sample is a good approximation 
of the one used in Ref. 4, all of the MARK I data including those measurements 
of the charm-threshold region are included and the Crystal Ball data and LEP 
measurement of ^^(MJ) are excluded. The resulting value of Att^°j*, 0.0231, is 
close to the value given in Ref. 4 of 0.0233±0.0009. 

The updated result given in Ref. 5 was derived by excluding the high energy 
MARK I data and including the Crystal Ball measurements. We simulate this re- 
sult by substituting the Crystal Ball data for the 44 high energy MARK I measure- 
ments in the equally-weighted analysis (the charm threshold data are still included 
and the LEP measurement of (^^(M^) is still excluded). The resulting value of 
Aa'j^^ , 0.0226, is close to the value given in the second publication of Ref. 5 of 
0.0228±0.0009. The best estimate of Rfit resulting from the equally-weighted fit 
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is shown as the dashed curve in Fig. 2. In addition to the difference in the charm 
threshold region, the equally- weighted fit prefers larger Rhad values in the region 
W = 1-3.6 GeV and in the region I^ > 10 GeV. In the W = 2.6-3.6 GeV region, 
the 15% normalization uncertainties of the DASP and PLUTO measurements pull 
the correctly-weighted fit to smaller Rhad values than those preferred by the 20% 
MARK I and 772 measurements. Since the 772 data extend to 1.42 GeV and 
are correlated by the large normalization uncertainty, the correctly-weighted fit 
is pulled to smaller Rhad values in the 1-2.6 GeV region. The use of the precise 
determination of Rhad sd, W = Mz is responsible for the Rfn difference in the 
high energy region. The Rf^t differences in the W regions 1-3.6 GeV, 3.6-6.0 GeV, 
and 8.0 GeV-00 lead to Aahad differences of 0.0003, 0.0004, and 0.0005, respec- 
tively. The contribution of the three ip resonances in the charm threshold region is 
discussed in Section 2.4 and is found to be 0.0002, approximately one half of the 
discrepancy in the 3.6-6.0 GeV region. We conclude that we have identified the 
largest part of the discrepancy and have validated our technique. 

The experimental uncertainty given in equation (20) is also smaller than the 
corresponding one given in Ref. 4. The authors of Ref. 4 evaluated their uncertainty 
by assigning normalization uncertainties in several W regions and adding them in 
quadrature. They assessed the following uncertainties: 20% in the region 1 GeV< 
W < 2.3 GeV, 10% in the region 2.3 GeV< W < 12 GeV, and 3% in the region 
W > 12 GeV. Note that the normalization uncertainty changes discontinuously at 
the boundaries of each region. The large uncertainties at small W are moderated 
by the fact that most of the integral accrues in the high energy region. Our use of 
RhadiMz) essentially eliminates the experimental uncertainty in the region W > 
15 GeV. The Crystal Ball data constrain the uncertainty in the region between 
c-quark threshold and b-quark threshold to be approximately one half of that 
assumed in Ref. 4. Below charm threshold, the normalization uncertainty increases 
from 15% near 3 GeV to 20% below 2 GeV. Since we require that Rj^t be a 
continuous (and reasonably smooth) function, the normalization uncertainty is 
required to vary smoothly as W decreases. The better measured larger W regions 
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constrain the uncertainty at the smaller W regions. 
2.3 The tt+tt" and K^K~ Final States 

The processes e+e" — >• tt+tt" and e+e" — )■ K^K~ are described by the elec- 
tromagnetic form factors, F^(s) and Fj{{s), which are related to the hadronic cross 
section ratio Rhad foi" each process as follows, 

Rllfis) = \\FAs)\'(il RLTis) = \\FK{s)\'fK. (21) 

where /3^ and l3j{ are the velocities of the final state particles in the e+e" center- 
of-mass frame. It is clear that measurements of the form factors are equivalent to 
measurements of Rhad- 

Measurements of the square of the pion form factor I^ttP have been per- 
formed by the OLYA,"^ CMD,"^ TOF,^'^ NA?,"^ /iTr,'"^ MEA,'"^ M2N,'"\ DMl,'"^ 

[I3l 

and DM2 Collaborations and are shown in Fig. 3. The error bars include the 
normalization uncertainties which range from about 2% in the region around the 
(dominant) p resonance to about 12% at IT ~ 2 GeV. 

The data are fit to a function which is a sum of the Gounaris-Sakurai form 
used by Kinoshita, Nizic, and Okamoto and three resonances, 

"" ^^ Al + A2q^ + fis) ^Z^^s-ml + imj^rr ' ^ ^ 

where: A\ and A2 are free parameters; m^ is the pion mass; q and /(s) are defined 
as follows. 



s/4 — ml 



fis) ^ - 



1 r o s 



m 



TT L '^ 3, 



2q 



3 



+ — ^In 



V^ + 2^' 



2m 



TT 






3 (23) 



and where m„, F„, B^, and C„ are the mass, width, amplitude, and phase of 
each resonance. The mass and width of the first resonance were set to those of 
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the uj(782). All other parameters (12 in total) were allowed to vary. The weight 
matrix of the fit was constructed by assuming that all normalization uncertainties 
are 100% correlated (see equation (13)). The result of the fit is shown as a solid 
line in Fig. 3. The fit preferred a resonance of width 0.36 GeV at mass 1.2 GeV 
and a second resonance of width 0.16 GeV at mass 1.7 GeV. The fit quality is 
found to be excellent (xV^of = 116.7/128). 
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Figure 3. Measurements of [.F^Cfy)!^ are compared with the best fit which is shown 
as a solid line. The error bars include normalization uncertainties. 

To evaluate the sensitivity of the result to the parameterization, the complete 
function used by the authors of Ref. 45 was also fit to the data. This function did 
not fit the newest data from DM2 (at large W) as well as our chosen form. Both 
functions were used to evaluate equation (5) from s = 4m^ to s = 4 GeV^ (where 
|F^P is measured to be very small). We find the Tr"*"?:" contribution to Aa^adiM'^) 
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to be 

Aafj'iM'^) = 0.003087 ±0. 000051 (exp) ± 0.000121 (param). (24) 

The two techniques for the estimation of the experimental uncertainty (discussed 
in Section 2.2) yielded consistent results in this case. 

The result given in equation (24) differs from the corresponding result given 
in Ref. 4 mostly because of the inclusion of the large-W^ DM2 data which decrease 
more sharply with energy than the extrapolated tail of the function that was fit to 
the lower-energy data. The uncertainty quoted by the authors of Ref. 4 corresponds 
to the fractional uncertainty on the leptonic width of the p. 

Measurements of the square of the kaon form factor \Fj{\'^ have been per- 
formed by the OLYA,"'^ CMD,"'^ MEA,'"^ DMl|"^ and DM2^"^ Collaborations 
and are shown in Fig. 4. The data span the (j){1020) resonance and continue to 
W = 1.8 GeV where R^^/^ is less than 0.01. The normalization uncertainty 
on the CMD measurements is 6%. The other groups do not report normalization 
uncertainties. Early \F^\ measurements suffered from the same problem of un- 
reported normalization uncertainties. A bit of historical research shows that the 
normalization uncertainties were usually not included in the point-to-point errors. 
We therefore arbitrarily assign a 20% systematic normalization uncertainty to all 
unreported cases. The data and total uncertainties are shown in Eig. 4. 

The data are fit to a function which is a sum of a Breit-Wigner resonance with 
an energy-dependent width for the (p and four resonances, 



A, ^ S„e^^» 



Fj,{s) = ""'. ^ , , +y " . ^ , (25) 



where: Ai is the amplitude of the cp; rn^ is the mass of the (/)(1020); m„, E„, B^, 
and Cn are the mass, width, amplitude, and phase of the resonances. The energy- 
dependent width Ta{s) is assumed to consist of contributions from the K~^K~ , 
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KiKs, and Stt final states, 



r^^) 




497 J^+ifL , 347^^^ 



+ 0.156Gf^(s) 



(26) 



where: V^i is the nominal value of the (f) width, /9+(s) = Wl — Arn^j^^l s is the 
velocity of the charged kaon, [3q[s) = -\/1 ~ 4m?^o/s is the velocity of the neutral 

- rn^i and is 



-i4> 



kaon, and G'^^{s) is a function which is normalized to unity at 



proportional to the decay rate for </)—)■ 3 tt assuming pvr dominance. 



[46] 
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Figure 4. Measurements of |_Fft-(iy)p are compared with the best fit which is 
shown as a solid line. The error bars include normalization uncertainties. 

The masses and widths of the first two resonances were set to those of the 
p(770) and tj(782). Following the procedure of Ref. 18, the amplitude ratios BilAi 
and B2/A1 were constrained to the measured values and the phases were set to 
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zero. The mass, width, and ampHtude of the (j) were allowed to vary. The masses, 
widths, amplitudes, and phases of two larger mass resonances were free parameters. 
The weight matrix of the fit was constructed by assuming that all normalization 
uncertainties are 100% correlated (see equation (13)). The result of the fit is shown 
as a solid line in Fig. 4. The fit preferred a resonance of width 0.15 GeV at mass 
1.39 GeV and a second resonance of width 0.22 GeV at mass 1.65 GeV. The fit 
quality is found to be adequate (% /dof = 74.1/49). 

To evaluate the sensitivity of the result to the parameterization, a second fit 
was performed with the amplitudes and phases of the p and co allowed to vary as 
free parameters. The fit quality improved marginally (x^/dof = 69.6/45). Both 
functions were used to evaluate equation (5) from s = 4:m^^+ to s = 3.24 GeV^. 
We find the K^K~ contribution to Aa^adiM'^) to be 

Aa£7''"(M|) = 0.000311 ± 0.000030(exp) ± 0.000009(param). (27) 

The two techniques for the estimation of the experimental uncertainty (discussed 
in Section 2.2) also yielded consistent results in this case. To account for the poor 
fit quality, the experimental uncertainty has been scaled by the factor 1.23. 

2.4 The Resonances 

The resonances comprise the remaining portion of the total e+e" cross section. 
The total cross section for each resonance can be represented by a relativistic Breit- 
Wigner form with energy-dependent widths, 

_ 127r sT,,{s)Ttot{s) 

where: m. Tee, and Ttot ai'e the mass, electronic width, and total width of the 
resonance. In order to incorporate the Breit-Wigner cross sections described by 
equation (28) into equation (5), it must be scaled to the electromagnetic point 
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cross section, o^^{s) = 4:7Ta {s)/3s, yielding the following expression, 



oo 



.res . 2^ _ ^p f , (^resis) 

which has the slightly unpleasant feature that it incorporates a{s), the quan- 
tity that we are attempting to evaluate, into the integrand. Equation (29) is 
often written in the approximation, a{s) ~ ao- Unfortunately, this overestimates 
AahadiMz) by: 3% at the cj(783), 5% at the J/V'(15), and 7% at the T(15). 
To avoid this problem, we use the Aaf^adis) parameterization given in Ref. 4 to 
generate a first-order estimate of a{s) for use in equation (29). 

Equation (29) is evaluated for the tj(782), 0(1020), '0-family, and T-family 
resonances by performing a Simpson's rule integration over the interval m — GOT tot 
to m -|- 60Ttot (the lower limit of the co integration is the threshold for Stt decay). 
The correction for electromagnetic decays is performed using the best estimate of 
the function Rfit(s) determined in Section 2.3. The energy-dependent electronic 
widths and the hadronic widths of the ip and T resonances are assumed to scale 
as ^/s, 



r{s) = ^Eo, (30) 

m 

where m is the mass of the resonance and Eq is the nominal value of the width. The 
energy-dependent total hadronic width of the <f){1020) is given by equation (26). 
The width E^^^ for the (j) is adjusted to exclude the K'^K~ final state (discussed 
in Section 2.4). The energy-dependent total hadronic width of the 6u(782) is given 
by the following expression which assumes that all final states are tt+tt", 7r*^7, or 
7r"'"7r~7r , 



r^(s) = r 



' ^'S 



m 



u> 



... rM , .nor (1-<A)^ 



0.022 'J^ ', +0.085 



/5^("^S) ' (i-"^^/"^S)^ 



+ 0.893G^,(s) , (31) 



where: ni^ is the mass of the tj, E^, is the nominal value of the uj width, /3,r('5) = 
yl — \m^ls is the velocity of the charged pion, and Gg^(s) is a function which is 
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normalized to unity at s = m^ and is proportional to the decay rate for u ^ Sty 
assuming a constant matrix element (phase space weighting). 

The masses and widths used to evaluate equation (29) are taken from the 1994 
Review of Particle Properties . The results are listed in Table 1 along with those 
derived in Sections 2.3 and 2.4. The experimental uncertainties are evaluated by 
assuming that the uncertainties on the masses, total widths, electronic widths, and 
relevant branching ratios are uncorrelated. The parameterization uncertainties are 
evaluated by repeating the calculation with a constant-width, constant-mass Breit- 
Wigner cross section. 

Table 1: Summary of the various contributions to Aaf^^d- 



Contribution 


W Region (GeV) 


A«w(M|) 


^{^Oihad)exp 


H^'^had)param 


Continuum 


l.O-oo 


0.021428 


0.000724 


0.000150 


7r+7r~ 


0.280-2.0 


0.003087 


0.000051 


0.000121 


K+K- 


0.987-1.8 


0.000307 


0.000025 


0.000009 


Resonances 


Ja) 


0.000305 


0.000010 


0.000003 


fi 


^(b) 


0.000304 


0.000011 


0.000004 


II 


ip (6 states) 


0.001106 


0.000059 


0.000023 


II 


T (6 states) 


0.000118 


0.000005 


0.000003 


Total 




0.02666 


0.00072 


0.00019 



(''^Doesn't include tt+tt final states. 
^^ Doesn't include K'^ K~ final states. 

The only resonance entry in Table 1 that is directly comparable to a corre- 
sponding result in Ref. 4 is T-family result which agrees well despite the use of 
a^{s) in equation (29). The sum of our K'^ K~ and (j) entries is larger than the 
corresponding (f) result of Ref. 4 by 16%. Taking the correction for a{s) into ac- 
count, this implies that continuum K'^K~ final states contribute about 20% of the 
(j) contribution to /S.aj^a.di^'z) ■ ^^^ ip-i-cLrai\y result is larger than the result given 
in Ref. 4 by 2% even though it includes three additional states. As a cross check. 
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a repetition of the calculation with T'^^^ = Ttot, 0((s) = otQ, the constant-width 
and mass Breit-Wigner function, and 1988 values of the resonance parameters did 
agree very well with the numbers given in Ref. 4. 

2.5 Final Result 

The various contributions to Aa/jaj(Mj) are summarized and summed in Ta- 
ble 1. The resulting value, 

Aahad{M^) = 0.02666 ± 0.00075, (32) 

differs by 2.6 (Ref. 4) standard deviations from the result given in Ref. 4 and by 1.9 
standard deviations from the updated result given in Ref. 5. Including the leptonic 
contribution, we find a~ (MJ) to be, 

a-i(M|) = 129.08 ±0.10, (33) 

where the uncertainties on the lepton masses contribute negligibly to the total 
uncertainty. 

3. Interpretation 

Since most electroweak calculations are based upon input parameters (tto, Mz, 
and the Fermi coupling constant, Gp) that are measured at very different scales, 
the quantity a{M'^) enters into the calculation of most electroweak observables. To 
understand the effect of our result upon the Standard Model predictions for various 
electroweak observables, we have used the value of a(Mj) given in equation (33) 
with the ZFITTFR 4.8 program of Bardin, et al. to calculate: the mass of 
the W boson (M^), the width of the Z boson (F^), the ratio of the hadronic 
decay width of the Z to the (single species) leptonic decay width (i?/), the tree- 
level total hadronic cross section at the Z pole (crj|^j), the effective weak mixing 
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angle at the Z pole (sin ^p|;), the ratio of the hh decay width of the Z to the 
hadronic width (i?^), the ratio of the neutral and charged current cross sections in 
neutrino nucleon scattering {Rv\ and the weak charge of the Cesium nucleus as 
measured in atomic parity violation experiments (Qvf(Cs)). The default value for 
Aahad{Ml) in ZFITTER is the result of Ref. 4 with the result of Ref. 5 available 
as an option. The shifts in the observables at a top quark mass (mt) of 175 GeV 
and a Higgs boson mass (mj/) of 300 GeV from those calculated with the default 
value of Att/jaj(M|) are listed in Table 2. The shifts are also normalized to the 
current world-average experimental uncertainties on each quantity. It is clear that 
the interpretation of the current measurements of sin 6^ is most affected and that 
the interpretation of the Tz and Ri measurements is moderately affected. 

Table 2: Shifts in the predicted values of various electroweak observ- 
ables using the value of Aa^adiM'^) given in equation (33) with the 
ZFITTER 4.8 program of Bardin et al. The shifts are calculated at 
mt = 175 GeV and nin = 300 GeV. 



Observable 


Shift wrt Ref. 4 [Ref. 5] 


Shift/Exptl Error 


Mw 


+40 [+30] MeV 


+0.22 [+0.17] 


Tz 


+ 1.9 [+1.4] MeV 


+0.50 [+0.37] 


Ri 


+0.013 [+0.010] 


+0.33 [+0.25] 


<d(Z) 


-0.004 [-0.003] nb 


-0.04 [-0.03] 


sin2 ^^ 


-0.00074 [-0.00055] 


-1.8 [-1.4] 


Rh 


-0.00003 [-0.00002] 


-0.015 [-0.010] 


R^{nc/cc) 


+0.0005 [+0.0004] 


+0.18 


QwiCs) 


+0.163 [+0.122] 


+0.09 [+0.07] 



The effect on the interpretation of sin 6^ can be made more clear by calcu- 
lating the allowed range for the currently favored value of the top quark mass, 
nit = 174: GeV, as the Higgs boson mass varies from 60 GeV to 1 TeV, 

0.2306, niH = 60 GeV 
sin^^^ = { 0.2315, mn = 300 GeV (34) 

0.2322, niH = 1 TeV. 
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The current determination of sin 6^ by the LEP Collaborations, 0.2321zb0.0004, 
is consistent with the CDF top mass value and a heavy Higgs boson. The value of 
sin 6^ extracted from the left-right cross section asymmetry in Z production by 
the SLD Collaboration, 0.2294±0.0010, is smaller than the light Higgs value by 
1.2 standard deviations. 



Table 3: The results of global MSM and S-T-U fits to the electroweak 
observables listed in Table 2. 



Parameter 


a-HMl) 




128.80 ±0.12 


128.87 ±0.012 


129.08 ±0.10 




Standai 


■d Model Fit 




mt (GeV) 


179 fi+10.2+16.3 
J^ '^■'^-10.8-18.1 


if^q 7+10.3+16.4 
J^"^-' -10.9-18.3 


Ifil 4+10.4+17.0 


as 


0.125±0.005 


0.125±0.005 


0.125±0.005 


xVdof 


22.3/12 


21.6/12 


20.0/12 




S-' 


T-U Fit 




S 


-0.37±0.23 


-0.31±0.23 


-0.16±0.22 


T 


+0.27±0.23 


+0.27±0.23 


±0.28±0.23 


U 


-0.20±0.56 


-0.18±0.56 


-0.14±0.56 


as 


0.124±0.005 


0.124±0.005 


0.123±0.005 


xVdof 


11.4/9 


11.5/9 


11.6/9 



We have performed a fit of the ZFITTER model to the current best measure- 
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ments of all of the quantities listed in Table 2 (the seven different observables that 
determine sin 6^ are entered separately). The fit was performed with the parame- 
ter a~^{M'^) constrained to 128.80±0.12 (corresponds to the value of Aa^adiM'^) 
given in Ref. 4), 128.87±0.12 (corresponds to the value of Att7jad(Mj) given in 
Ref. 5), and 129.08±0.10. The full correlation matrix for the LEP measurements 
is included in the fit. The top quark mass and the strong coupling constant were 
allowed to vary as free parameters in a series of three fits with the Higgs boson 
mass set to: 60 GeV, 300 GeV, and 1 TeV. The results are listed in Table 3 for 
the rriH = 300 GeV case. The extracted value of mt is sensitive to the choice of 
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Higgs mass. The range of sensitivity is indicated by the second set of errors as mjj 
is varied from 60 GeV (lower value) to 1 TeV (upper value). Note that the use of 
our new value of a{M'^) decreases the extracted value of m^ by 11 GeV [8 GeV] 
as compared with the result of Ref. 4 [Ref. 5] and slightly improves the quality of 
the fit. 

We have also performed a fit of the Peskin-Takeuchi S-T-U parameterization 
to the best current measurements of the quantities listed in Table 2. The param- 
eters 5, T, [/, and ttjj(Mj) were allowed to vary as free parameters with a~^{M'^) 
constrained to 128.80±0.12, 128.87±0.12, and 129.08±0.10. The results are listed 
in Table 3 for the reference values of mt and rriH taken to be 150 GeV and 1 TeV, 
respectively. The use of the new value of a{M'^) shifts the extracted value of S by 
+0.21 [+0.16] as compared with the result of Ref. 4 [Ref. 5]. 

4. Conclusions 

We have reevaluated the hadronic part of the electromagnetic vacuum expecta- 
tion value using the standard dispersion integral approach that utilizes the hadronic 
cross section measured in e+e" experiments as input. Previous analyses are based 
upon point-by-point trapezoidal integration which has the effect of weighting all 
inputs equally. We use a technique that weights the experimental inputs by their 
stated uncertainties, includes correlations, and incorporates some refinements. We 
find the hadronic contribution to Aa{M'^) to be, 

AahadiMz) = 0.02666 ± 0.00075, 

which leads to a value of the electromagnetic coupling constant at s = M^, 

a-i(M|) = 129.08 + 0.10. 

Our value of a{M'^) shifts the predicted values of a number of electroweak 
observables. The most affected is the effective weak mixing angle at the Z pole 
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which shifts by —0.0007 [—0.0006] from that predicted using the Aa^adiMz) value 
of Ref. 4 [Ref. 5]. The best estimate of the top quark mass extracted from a global 
fit shifts by —11 GeV [—8 GeV] and the best estimate of the Peskin-Takeuchi S 
parameter shifts by +0.21 [+0.16]. 

We note that the current generation of sin 6^ measurements are likely to sat- 
urate the ±0.00026 uncertainty due to the ±0.10 uncertainty on a~^(M'^). The 
best hope for improvement is for the BES Collaboration to make improved mea- 
surements of Rhad ill fhs region W = 2 — 3.6 GeV. A modest set of measurements 
in this region with a normalization uncertainty SR^ad/Rhad ^ 5% would reduce 
the current uncertainty by a factor of two and would eliminate it as a limitation. 
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